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The complete analytical solutions of the geodesic equation of massive test particles in higher 
dimensional Schwarzschild, Schwarzschild-(anti)de Sitter, Reissner-Nordstrom and Reissner- 
Nordstr6m-(anti)de Sitter space-times are presented. Using the Jacobi inversion problem restricted 
to the theta divisor the explicit solution is given in terms of Kleinian sigma functions. The derived 
orbits depend on the structure of the roots of the characteristic polynomials which depend on the 
particle's energy and angular momentum, on the mass and the charge of the gravitational source, 
and the cosmological constant. We discuss the general structure of the orbits and show that due 
to the specific dimension-independent form of the angular momentum and the cosmological force 
OO ■ a rich variety of orbits can emerge only in four and five dimensions. We present explicit analytical 

' solutions for orbits up to 11 dimensions. A particular feature of Reissner-Nordstrom space-times is 

' that bound and escape orbits traverse through different universes. 
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I. INTRODUCTION AND MOTIVATION 

O ! The idea of higher dimensions has a long history. Already in the 1920s Kaluza and Klein unified electromagnetism 
and gravity by introducing a 5th dimension, while today string theory, as a promising candidate for the quantum 
theory of gravity and for the unification of all interactions, requires higher dimensions for its mathematical consistency. 
Following Kaluza and Klein the higher dimensions may be considered small and compact and therefore not observable. 
Alternatively, in braneworld scenarios, large higher dimensions are invoked, where only gravity is allowed to penetrate 
, the extra dimensions [l|, Q , while the other fundamental interactions are considered to be trapped on the 3-brane 

^ ' forming our universe. 

OO , In such braneworld models the four-dimensional gravitational coupling G, corresponding to the Planck scale of 
10^^ GeV, does not set the fundamental scale. Instead, the fundamental scale is assumed to be on the order of the 
[ electroweak scale, and the extra dimensions then serve to solve the hierarchy problem by invoking either flat compact 
. extra dimensions as in the Arkani-Hamed-Dimopoulos-Dvali model [l| or warped extra dimensions as in Randall- 
■ Sundrum models ■ Both type of models predict the production and evaporation of mini black holes in high-energy 
collisions Q at energies in the TeV range accessible at the LHC at CERN or by observing high-energy cosmic rays 
in the atmosphere [J Q . 

The study of higher dimensional space-times is also of relevance from another point of view. The AdS/CFT 
correspondence Q suggests a duality between a string theory involving gravity in Anti-de Sitter (AdS) space, and 
a conformal field theory (CFT) living on its boundary. Study of the phases of higher dimensional AdS black holes, 
• for instance, can then shed light on the properties of strongly coupled gauge theories. 
d One of the first papers dedicated to the problem of finding solutions of the Einstein equations in higher dimensions 

both in vacuum and coupled to matter fields was by Tangherlini, who presented the Schwarzschild and Reissner- 
Nordstrom solutions in n dimensions 3J- Two decades later Myers and Perry generalized the rotating Kerr black 
hole solutions to higher dimensions. The metric of a rotating black hole with a cosmological constant was derived 
in 5 dimensions and in TTl for arbitrary dimensions, while the general Kerr-NUT-(A)dS solution was presented 
. A review of black hole solutions in higher dimensional vacuum gravity was given by Emparan and Reall [Tsj . In 
contrast, the higher dimensional charged rotating black holes of Einstein-Maxwell theory have not yet been obtained 
in closed form [14|. Perturbative Einstein-Maxwell rotating black holes were studied in [H, [H, [13, [3, while non- 
perturbative black hole solutions were obtained numerically in [l^ and a finite cosmological constant was included 
in [21], (see [11] for a recent review) . 
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To understand the physical properties of solutions of the gravitational field equations it is essential to study the orbits 
of test particles and light rays in these space-times. On the one hand, this is important from an observational point 
of view, since only matter and light are observed and, thus, can give insight into the physics of a given gravitational 
field [24]. The study of the motion of test particles in gravitational fields also has significant practical applications. 
On the other hand, this study is also important from a fundamental point of view, since the motion of matter and 
light can be used to classify a given space-time, to decode its structure and to highlight its characteristics. 

When integrating geodesic equations two fundamental questions arise. These concern the separability of the 
Hamilton-Jacobi equation and the possibility to solve the equations of motion analytically. The first question is 
related to the symmetries of a given space-time. Continuous isometries preserving a metric in curved space-times 
are generated by Killing vector fields. For each Killing vector field there exists a conserved quantity. For example, 
when all metric components do not depend on time and the azimuthal angle, energy and angular momentum are 
conserved during a particle's motion. In 1968 Carter (26j discovered a new separation constant in the Kerr space-time 
and, thus, proved the separability of Hamilton-Jacobi equations. This peculiarity of the Kerr metric is associated 
with hidden symmetries associated with the Killing- Yano [13] tensor. Geodesic equations are completely integrable, 
and Hamilton-Jacobi, Klein-Gordon and Dirac equations are separable in space-times with a metric admitting a 
conformal Killing- Yano (CKY) tensor. A CKY tensor generates a separation constant proportional to a quadratic 
an gula r momentum of a particle. Hidden symmetries and the separation of variables in higher dimensions are discussed 
in 

The existence of analytical solutions for any kind of problem is not just an academic question. In fact, analytical so- 
lutions offer systematic applications as well as a frame for tests of the accuracy and reliability of numerical integrations. 
In 1931 Hagihara [s^ first analytically integrated the geodesic equation of a test particle motion in Schwarzschild 
gravitational field. This solution is given in terms of the Weierstrass p-function. With the same mathematical tools 
one can solve the geodesic equation in a Reissner-Nordstrom spacetime [33|. The solutions of the geodesic equation 
in a Kerr and Kerr-Newman space-time have also been given analytically (see [ssj for a survey). Recently two of 
us found the com plet e analytical solution of the geodesic equation in Schwarzschild-(anti) de Sitter space-times in 
4 dimensions (3^. l35l|. These calculations are based on the inversion problem of the hyperelliptic Abelian integrals 
[S^l- The equations of motion could be explicitly solved |33, 11^ by restricting the problem to the theta divisor. This 
procedure makes it possible to obtain a one-parameter solution of the inversion problem. The approach was suggested 
by Enolskii, Pronine and Richter who applied this method to the problem of the double pendulum [38| . 

Our interest here focuses on the motion of test particles in higher dimensional spherically symmetric space-times. 
In Table [J we characterize space-times with respect to the possiblity to solve the geodesic equation analytically. 
Space-times whose geodesic equations can be integrated by elliptic functions are marked by a +. The underlying 



Space-time — — — ^ ' 
^ — — " Dimension 


4 


5 


6 


7 8 9 10 11 d > 12 


Schwarzschild 


+ 


+ 


e 


+ - e - e - 


Schwarzschild-de Sitter 


e 


+ 




+ - © - ffi 


Reissner-Nordstrom 


+ 


+ 




ffi 


Reissner-Nordstrom-de Sitter 


® 


+ 




ffi 



TABLE I: Characterization of higher dimensional spherically symmetric space-times with respect to the integrability of the 
geodesic equation. The + denotes the solvability by elliptic functions, ffi the solvability in terms of hyperelliptic functions, and 
a — indicates that it is not analytically solvable by any known method. 

polynomials P{r) in these equations can be reduced by a simple substitution r = or r = 1/v^ to polynomials of 
degree 3 or 4. Polynomials of 4th order can be subjected to a subsequent substitution x = 1/u + r4 which reduces 
Pi{x) to a cubic expression (here is a zero of Pi{x)). The complete integration of the geodesic equation with an 
underlying polynomial of order 3 can be performed in terms of elliptic functions. As an example consider the geodesic 
equation in the Schwarzschild space-time in 7 dimensions 



^y^MJ ^here Pe{r) ^ (E^ - 1)7-' - r^L' + + L^l , (1) 
dip J r'^L^ 

with the conserved quantities E and L and rg — 2M 42]. Introducing the new variable x = r'^ reduces equation ([T]) 
to a new one with a cubic polynomial Paix) 

(7t \ ^ 4 

_ ^—P^ix) where P^ix) ^ (E^ ~ l)x'' - x"" + 4x + L^l . (2) 
dip J 
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The solutions of equation ^ are given in terms of the Weierstrass p-function. 

The symbol ® indicates space-times whose geodesic equations contain polynomials of order 5 (or polynomials 
reducible to a polynomial of the 5th degree). These equations are hyperelliptic and can be solved by the method 
suggested in [3J, |3a] . In these papers the method was applied to the integration of geodesies in the Schwarzschild-de 
Sitter space-time in 4 dimensions and in the Schwarzschild space-time in 6 dimensions. Here we study the remaining 
seven cases. The "— " indicates space-times whose equations of motion cannot be integrated analytically by any 
known methods, since the underlying polynomials are of 7th order or higher. 

We begin with a short introduction to the new method of hyperelliptic functions used for the integration of the 
equations of motion and apply it first to the Schwarzschild space-time in 9 and 11 dimensions in Section [TVl Then we 
continue with the Schwarzschild-(anti-)de Sitter space-time in both 9 and 11 dimensions in Section |Vl the Reissner- 
Nordstrom (RN) space-time in 7 dimensions in Section IVII and the RN-(anti-)de Sitter space-time in 4 and 7 di- 
mensions in Section IVIII We also include the solution of the geodesic equation in the 4-dimcnsional RN-(anti-)de 
Sitter space-time, which also requires the knowledge of the theory of hyperelliptic functions. The resulting orbits are 
classified in terms of the energy and the angular momentum of a test particle, the charge of the gravitational source 
and the cosmological constant. 



II. THE GEODESIC EQUATION 



A. The space-times 



We restrict ourselves to static spherically symmetric vacuum solutions of the Einstein equations in an arbitrary 
number of dimensions d 



with 



ds^ = f{r)de - f{r)-Hr^ - r^dnl_^ , (3) 



where M and q are the mass and charge, respectively, of a gravitational source, A is a cosmological constant, rg — 2M , 
and dVL\ — dip^ and dV,^_^_^ — dOf + sin^ 9idflf for i > 1. These solutions are uniquely characterized by their mass, 
charge and cosmological constant ^] . 



B. The geodesic equation 



The motion of a particle is given by the geodesic equation 

DyV = , (5) 

where v is the 4-velocity of the particle and D is the covariant derivative based on the space-time metric. We also 
use the normalization g{v,v) = e where e = 1 for massive particles and e = for light. 

Due to spherical symmetry of the space-time the motion of a test particle can be restricted to the equatorial 
plane defined by 9i = 7r/2, i > 1. There are two conserved quantities: the dimensionless energy E, and the angular 
momentum L with the dimension of length (both are normalised to the mass of the test particle) 

Equation ([5|) combines the equations describing the dynamics of the particle's motion 

E' - m ( e + !i] (8) 



ds J V 

dr_y _ fH^ 

dt ^ 
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We also introduce the effective potential Vcs which can be read off 



Kff = ( 1 - ( - 



d-3 



2Ar 



(d- 2) Vr 



g\2(<i-3) 



It is obvious that at any horizon, which is defined by /(r) = 0, the effective potential vanishes. 

For convenience we rewrite ([7]) with new dimensionless parameters and a dimensionless coordinate f 



and obtain 



dip 



„d-3 



3 {d~l){d~2) f2(d-3) 



1 

Af2 



(10) 



(11) 



(12) 



The right hand side of this equation can be written in the form Pn{f)lf™' i where P„ is a polynomial of order n, see 
Table ni 





?7 = 




A = 
A / 


(dr_y _ Pd-1 
\ dip 1 — 5 

(dr.] - ^rf+l 


( dr\ P2(d-2) 
\d^ ) f2(d-i) 

( df \ ^2(d-l) 
\d^p ) f2{d-i) 



TABLE II: Order of the polynomials appearing in the equation of motion p2p in various d-dimensional space-times. 



Depending on the space^imc under consideration the equation of motion p2p after some transformation f = x{^) 
can be reduced to one of the two forms 



dx 
dip 



(13a) 
(13b) 



where ^5(2;) = ac^x^ + a^x'^ + a^iX^ + a2x'^ + aix + ag. 

Because of (jl2p space-times with odd number of dimensions give polynomials P„ containing even powers of f only. 
Therefore we simplify ((T2)) by introducing u — \ jr^ what reduces the degree of P„ by a factor of 2 



1 ( du\. 



4: \ dip J 



2A \ 2eAA 



+u X{fi - e) + 



(d-l)(d-2)/ (d-l)(d-2) 



(14) 



The number of positive real zeros of the polynomials on the right hand side of and (|12p coincide. Space-times 
with an even number of dimensions, as e.g. 4-dimcnsional Reissner-Nordstrom-de Sitter space-time considered in 
Section IVlI Al contain polynomials of odd degree so that this substitution is not useful. 



C. The topology of the orbits 

Owing to the square on the left-hand-side of ([S]) a test particle is only allowed to move in the region restricted 
by the condition > Vcff- Therefore, the motion is essentially characterized by the number of real positive zeros 
the polynomials Pnif) possess. The number of real positive roots of a polynomial Pn{f) can easily be determined by 
using Descartes' rule of signs. 

If Pn(f) possesses no real positive zero and Pn{f) > 0, then the particle is coming from infinity and moves directly 
to the origin. Such an orbit we call a terminating escape orbit. For Pn{f) < no motion is possible. One real positive 
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zero means either that the particle starts at a finite coordinate distance and ends at f = 0, what we cah a terminating 
bound orbit, or it moves on an escape orbit with a finite impact parameter. For two real positive zeros we may have 
two cases: (i) if the polynomial is positive between the two zeros we have a periodic bound orbit like a planetary orbit, 
(ii) if the polynomial is negative between the two zeros we have an escape orbit or a terminating bound orbit. The 
structure of the orbits for more than two positive zeros can be obtained analogously. The actual orbit depends on the 
chosen initial position. 

Below we discuss in more detail the structure of the orbits for the various dimensions and the allowed values of the 
energy, angular momentum and cosmological constant. As one result we see that for all considered space-times the 
richest structure of orbits can be found in four and five dimensions. The reason for that is that only in these lower 
dimensions various terms combine such that more possibilities arise to obtain different signs for coefficients of the 
polynomial. For higher dimensions all terms disentangle and the structure of the orbits becomes very simple and is 
the same for all dimensions. That means that it is the particular structure of the angular momentum part and the 
cosmological part both being independent of the space-time dimensions, which single out the lower dimensions. 



1. Schwarzschild space-times 



For the geodesies in Schwarzschild space-times in arbitrary dimensions we have 

{£)-^, (H, - ly-' - r'-' + Ar- + 1) , (15) 
which yields the types of orbits listed in Table WUl 

dimension polynomial parameters orbits 

d = 4 A(/i — l)f + — f + 1 ^ < 1 periodic bound, terminating bound 

H > 1 escape, terminating bound, terminating escape 
d = 5 — l)f* + (A — l)f^ + 1 ^<1, A<1 terminating bound 

pi < 1, A > 1 terminating bound 

pi > 1, A < 1 terminating bound, escape, terminating escape 
pi > 1, A > 1 terminating escape 
d > 6 A(pi - l)f''^-^ - f^^'^ + Af^ + 1 pi < 1 terminating bound 

pi > 1 terminating bound, escape, terminating escape 



TABLE III: Types of geodesies in Schwarzschild space-time 



Only in four dimensions there are periodic bound orbits. There are no stable configurations in higher dimensions. 



2. Schwarzs child- de Sitter space-times 



The geodesic equation in Schwarzschild-de Sitter space-time in higher dimensions 



dr^^ 



dip 



2AA 



(d-l)(d-2) 



pd+l 



+ A(pi-1) + 



2A 



{d-l){d-2) 



(16) 



contains an additional term Af'^"''^ as compared to the Schwarzschild case. For large f this term is dominant. For 
small negative A all escape orbits turn into periodic bound orbits, and for small positive A some periodic bound orbits 
(slightly below pi = 1) will turn into escape orbits. A huge positive A might surpass the angular momentum barrier. 
We skip the detailed discussion of all subcases and just remark that only for d = A and 5 different shapes of orbits 
occur due to the mixing of different terms. Such mixing no longer occurs for d > 6. 
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3. Reissner-Nordstrom space-times 

The geodesies in RN space-times 

do not possess any terminating escape orbit. The reason is the charge which leads to gravitational repulsion near the 
singularity at the origin. Accordingly, we have either periodic bound or escape orbits. A particular feature is that the 
4 dimensional RN space-time is the only space-time which allows 2 different periodic bound orbits (see Table WV\ . 
Both periodic bound orbits are characterized by the same period which is a consequence of the fact that the solution, 
owing to the order of the polynomial, is given in terms of the Weierstrass p function. 



dimension 


polynomial 




parameters 






orbits 


d = A 


A(/i - l)f* + Af^ - (1 + Ary)f2 + f - r) 




M < 1 






2 periodic bound 








1 < ^ 






periodic bound, escape 


d = 5 


A(m - 1)^='' + (A - l)f^ + (1 - Ar?2)f^ - rf 


11 < 


1, A < 1, A?7^ 


> 


1 


not allowed 






M < 


1, A < 1, Xrf 


< 


1 


periodic bound 






M < 


1, A > 1, Xrf 


> 


1 


periodic bound 






M < 


1, A > 1, Xrf 


< 


1 


periodic bound 






M > 


1, A < 1, Xrf 


> 


1 


escape 






M > 


1, A < 1, Xrf 


< 


1 


periodic bound, escape 






M > 


1, A > 1, Xrf 


> 


1 


escape 






M > 


1, A > 1, Xrf 


< 


1 


escape 


d > 6 


A(^ - l)f2(''-2) - f2(<*-3) + Af'^-i 




/i < 1 






periodic bound 




+f'*-^ - Ar,^-3f2 - r^^-^ 




1 < ^ 






periodic bound, escape 




TABLE IV: Types of geodesies in Reissner-Nordstrom space-time 



^. Reissner~Nordstrdm-de Sitter space-times 
For the addition of a cosmological constant in the geodesic equation in a RN space-time 



d(p 



= Xf 



2Af 



?2 



„d-3 



3 (d-l)(d-2) f2(d-3) 



1 

Af2 



(18) 



the remarks made after p6p also apply. A point to add is that the periods of the two periodic bound orbits appearing 
in RN space-time now become different. 



D. Integration 



Separation of variables in p3)) yields the hyperelliptic integral 



or 



x'dx' 



dx' 



(19a) 



(19b) 



for the physical angle ip with xj^ — x{ipia) (the pair (a;in, y'in) is an initial point of a test particle motion). This approach 
allows us to obtain an explicit solution of the inversion problem by restricting the problem to the theta-divisor (zeros 
of the theta function). As a result we obtain the complete set of analytically given orbits of test particles. 
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There are two important points: The first is the ambiguous definition of the integrand because of two branches of 
the square root. The second is the periodicity of x{ip) following from the requirement of x{ip) to be independent of the 
integration path. That means that for a closed integration path 7 and tu = §^ ^p^( ') ^ ~ ^1 y'p^{ >) ' ''^spectively, 
the function x{ip) has to fulfill 

r x'dx' r dx' 

and, thus, has to be a periodic function with period to. 

Owing to some modifications compared to [35l | we shortly review the procedure of explicit integration. 



III. SOLVING HYPERELLIPTIC INTEGRALS 



We consider (fT9|) on the Riemannian surface — P5{x) of genus g — 2 and introduce a basis of canonical 
holomorphic and mcromorphic differentials dzi and dui, respectively, 

dx xdx , , 

dzi = — , , dz2 = — , , (21) 

aax + 2042;^ + Sa^x^ x'^dx 
dui = , dx , du2 — , , (22) 

and real 2uJij , 2?}ij and imaginary 2uj'^j , 2ri'^j period matrices 

2u}tj — j) dzi , 2uo[^ = (f) dzi , (23) 

J aj 

2riij — (p du. 



bi 



2?7^ = / du, (24) 



with i,j = 1, 2 and the canonical basis of cycles (ai, &i, 62)- These 4 closed paths defined on a Riemannian surface 
correspond to periods of a function defined on it (in our case, to the solution of (|19p V We also need the normalized 
holomorphic differentials 

dv={2u})-^dz, dz=['^''A. (25) 

\dz2 J 

The period matrix of these differentials is given by (I2, t), where r is a Riemann matrix defined by t = lu^^lu' . It is 
this T which contains all the information about the polynomial P5. 

The solutions of can be inferred from the solutions of the Jacobi inversion problem [s^] which solves 



"Pi 



dx dx 



xo y/Psixj Jxa \/Pb[x) 



, xdx f^^ xdx 
V2 = / . + / . (26) 



for X = (xi, in terms oi ip = {(fi, ^2)^ ■ This solution is given by 



4 

X\ + X2 = P22KV) I 

4 (27) 

XiX2^ Pl2{<?) , 

as 



where 



,^ d d , ,. (Ti{z)aAz) — a{z)aij(z) 
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are the generalized Weierstrass functions ((7^ = da/dzi), where 

cj[z)^Ce~^^'^^''^^({2Lo)-^z + Tg + h-T) , (29) 
(the constant C can be given expUcitly, see [i^) is the Kleinian a function based on the theta function 

It is the i9-function where the matrix r conies in. Note that the problem (|26p contains the problem (jl9p for X2 — » a;o 
and a restriction such that one of the two equations ((26| can be omitted. 

Without restriction we can choose xq = oo (this only induces a reformulation of Lp) and consider the limit X2 oo 
to obtain 

XiX2 Cr(<Poo)0-12(<^oo) - tTl((^oo)cr2(<^oo) .„,^ 

XI = hm ■ = 27^— — — -— , 31) 

X2^o:>Xi+X2 Cr|((^oo) - a-(c/?oo)0-22(¥'oo) 

where (/Joo — liniajj^oo dz ~ f^^ dz — dz, where x-m is an initial point of a test particle motion. 

The next step is to restrict Jacobi's inversion problem to a theta-divisor Oa'^o i the one-dimensional set of zeros of 

the theta function -i? ^(2a;)^^z + i^^r,-,; , where K^o is the Riemann vector associated with the base point xq [40| . 

One can check that (2Lu)^^(pao is an element of the theta-divisor 8^ with Kao = (1/2) + (1/2) Si and 
a{(poc) = 0. This allows us to write pip in a simpler form and to express one of the components of ipoo in terms of 
the other component. 

Now we can write the analytic solution of the geodesic equation as 

f (^) = xi^if )) = Xixi) - X . (32) 

Eqn. ([5^ is an explicit analytical solution of the geodesic equation (|12p which we will explicitly work out for the 
Schwarzschild, Schwarzschild-de Sitter, Reissner-Nordstrom and Reissner-Nordstrom-de Sitter space-times in higher 
dimensions. 

It is important to note that, depending on the type of the problem ()19p . there are two possible physical angles (p: 

= dz2 + (fiin , (33a) 
fzi = dzi + (pin . (33b) 

J 

This means that depending on the problem ^Pao, which is an element of the theta-divisor and, therefore, one- 
dimensional, will take a corresponding form 

Cr^^ dzi — f°° dzi\ , , 



2 , 



where (^J„ = fin + /^°° dz2 and (p[^ = ip\n + dz\. Note that in (foo.z^ the first component is a fimction of the 
second component and in V'oo.zi the second component is a function of the first component. 

Now we are going to explicitly work out the solutions for various higher dimensional space-times. 



IV. GEODESICS IN HIGHER DIMENSIONAL SCHWARZSCHILD SPACE-TIMES 



In this section we present and discuss the possible types of orbits in the space-time of a gravitating source described 
by the higher dimensional Schwarzschild metric given by ([3]) for A = and q — 0. We first address the effective 
potential V^s in these space-times and describe the classification of zeros of the polynomial P„ (f) in (|12p in terms of 
the parameters and A. We then present the explicit analytical solution of the geodesic equation for e — 1. 



(a) 



(A, /x)-plot 



(b) 



Potential for A = 0.15 



FIG. 1: 9-dimensional Schwarzschild space-time, (a) The zeros of P5 in a (A, /i)-diagram (tlie number of zeros is encoded in 



the gray scales: white = 0, light gray — 1, gray — 2 zeros), (b) The effective potential. For the corresponding orbits see Fig. [2] 



A. Schwarzschild space— time in 9 dimensions 

The equation of motion p4p for Schwarzschild space-time in 9 dimensions reduces to 

(^) =4u(u4 + Au3-u + A(//-l)) =4P5(m) • (35) 



The A — /i diagram in Fig. 1(a) shows the number of positive zeros of the polynomial P5{u) (or, equivalently, Pg(f)), 
depending on the parameters A and /i. Here the gray scale code is as follows: gray denotes two, light gray one, and 
white no positive real zero of P5{u). 



The effective potential pH)) is 



^cff=(l-^) (1 + ^) (36) 



With z = the extrema of the effective potential are given by = — 3Az — 4. The discriminant of this cubic 
equation is = 4 — A''. For D < there is always one positive real zero given by z = 2\/Acos 



, arccos — S= 

, V A-^ 



For £) = the positive zero is given by z = 4^/^ and Kfi (4^/^) = ||. For D > we have z = (2 + \/4 - A^) 

,1/3 



(2- V4^^) ' . For any value of A the second derivative of the potential at the corresponding value of z is negative 
which implies that the potential has a maximum. There are no stable circular or periodic bound orbits. 
Eq. ([55)1 is of type (|13bp and the solution ([5^ becomes 



'^(W = \ 7^ T , (37) 

with the angular coordinate ipz-^ of a test particle from (|33b[) and V'oo.zi from equation (|34bp : here is the starting 
point of the motion usually chosen to be a zero of Psix). 

The corresponding orbits of massive test particles are illustrated in Fig. [2] The chosen energies /i = 2.6 and /i — 3.9 



are shown as dashed lines in the effective potential Fig. 1(b) which swap out the regions where a particle is allowed 



to move. For both chosen energy parameters terminating bound and escape orbits exist. For /J, > 1 also escape orbits 
become possible. For higher energies both kinds of orbits become more deformed and approach the unstable circular 
orbit corresponding to the maximum of the potential. Then terminating escape orbits arise. 

Let ei . . . 65 be the roots of the polynomial of 5th degree Psiu). The deflection angle for an escape orbit of a test 
particle in Schwarzschild space-times in 9 dimensions is 

Ju^ v/4AR 1 -4^12 - ^ if 3fJe, > Ue 
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(c) fi = 3.9203: terminating bound orbit and escape orbit (d)/x = 3.92046627: terminating bound orbit and escape orbit 



5- 



(e) fj, = 4.0: terminating escape 
orbit 



FIG. 2: Orbits of a test particle in the 9-dimensional Schwarzschild space-time for A — 0.15. The black circle indicates the 
horizon. The plots are given in units of rs. 



where Ue is a zero of P5 (u) corresponding to an escape orbit and Sfte^ is the real part of a corresponding complex root 

a. 



B. Schwarzschild space— time in 11 dimensions 



For the Schwarzschild space-time in 11 dimensions the equation of motion p4p takes the form 

4m {u^ + Xu^-u + A(/i - 1)) = APeiu) . (39) 



The number of zeros of this Pg are shown in Fig. [3] in the (A, fi)-p\ot together with the effective potential 



1 \ / , 1 
Af2 



^cff = ( 1 - ) ( 1 + T-I ) (40) 



which indicates no qualitative difference between the Schwarzschild space-times in 9 and 11 dimensions. With z = f ^ 
the extrema of the effective potential are given by = — AXz — 5. The only possible real positive zero is given by 
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(a) 



(A, /x)-plot 



15 30 45 r 

(b) Potential for A = 0.4 



FIG. 3: 11-dimensional Schwarzschild space-time, (a) (A,^) plot (color code as in Fig. lamuSchwQD). (b) EfTective potential. 
For the corresponding orbits see Fig. |4l 
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(a) n = 2.1: terminating bound orbit and escape orbit (b) fi = 2.1760084521: terminating bound orbit and escape 

orbit 

FIG. 4: Orbits of a test particle in the Schwarzschild space-time in 11 dimensions for A = 0.4. In (b) the energy is close to the 
energy of the unstable circular orbit related to the maximum of the effective potential in Fig. |3][b). 



the expression: z ^ ^ i^^/k + i^-k + j j with ^ {27X^ + 3V81A4 + 375) ' - 5 {27X^ + SVSIA^ + 375) ' 
and describes a maximum. 

There are no stable circular or periodic bound orbits. A substitution u = ^ + uq, where uq is a root of Pe{u), 

transforms ((39)) to ^a:^^ — 'iP^{x) which is of type (|13ap . Then the solution of the geodesic equation for a test 
particle in an 11 dimensional Schwarzschild space-time is given by 



cri(Voo,^2) 



(41) 



with the angular coordinate ipz^ from (|33ap and ipoc,z2 from p4ap . is again an initial angle which can be chosen 
to be a zero of ^5(2;). The orbits for various /i are illustrated in Fig. |4]and exhibit the same qualitative behavior as 
the orbits in 9 dimensions. 

The deflection angle for an escape orbit of a test particle in the Schwarzschild space-time in 11 dimensions is 



cdx 



where Xe is a zero of ^5(2;) corresponding to an escape orbit and we took ug — here. 



(42) 
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FIG. 5: Potential for negative and positive cosmological parameter A. For the same value of a test particle energy /i it can be 
on a periodic bound orbit with perihelion and aphelion points A and B or on an escape one with a point of nearest approach 
A. 



V. GEODESICS IN HIGHER DIMENSIONAL SCHWARZSCHILD-(ANTI-)DE SITTER SPACE-TIMES 

We now consider timelike geodesies in Schwarzschild-(anti-)de Sitter (S(a)dS) space-time. A detailed discussion of 
geodesies in Schwarzschild-de Sitter space-time in 4 dimensions can be found in [35| . Here we restrict our consider- 
ations to 9 and 11 dimensions. The geodesic equations for 5 and 7 dimensions can be solved using the Weierstrass 
p-function. 



A. Schwarzschild— (anti-)de Sitter in 9 dimensions 

For d = 9 and 77 we obtain from (|14p the equation of motion in the S(a)dS space-time in 9 dimensions 

(1)^4 (.^ + A.^ - + (a(, - 1) + 4) ^ + = 4A (.) , (43) 

which is of type (jl3b[) . The solution of the geodesic equation in the 9 dimensional Schwarzschild-(anti-)de Sitter 
space-time is of the form ([57]) . The sign of the cosmological constant is crucial for the form of the orbits. As it can 
be seen from the effective potential 

shown in Fig. O for A < a particle necessarily follows a bound orbit, while for A > there are also escape orbits. 



Positive cosm 



constant 



For A > the positive zeros of P5{u) (and also of Pio(f)) are shown in Fig. 6(a) From this Figure as well as from 
the effective potential Vcs (|44p shown for A = 0.2 and A = 8.7 • 10^^ in Fig. 6(b) it is clear that only terminating 
bound, terminating escape and escape orbits are possible. The qualitative structure is mainly the same as for the 
pure Schwarzschild case except that for small energies the escape orbits show reflection at the cosmological constant 
barrier, Fig. |6(b)[ 



2. Negative cosmological constant 



For A < the positive zeros of P^^u) are shown in Fig. 8(a) For a negative cosmological constant the effective 
potential Vcff differs insofar as it tends to +00 for r — > 00. Therefore each orbit is a bound orbit, some are periodic, 



some are terminating. This is related to the dark gray region in Fig. 8(a) corresponding to three real positive zeros 



giving one terminating bound orbit and one periodic bound orbit, and the light gray region indicating one real positive 
zero leading to terminating orbits (see Fig. [5]). 
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(a) (A,Ai)-plot, A = 8.7 • 10- 




/i = 3.5 
fj. = 3.13 

Ai = 1.8 
Ai = 1.008 
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(b) Potential for A = 0.2, 
A = 8.7 ■ 10-5 



FIG. 6: 9-dimensional SdS space-time: (a) (A,/i)-plot. (b) Effective potential. For the corresponding orbits see Fig. [7] 



For a periodic bound orbit restricted to the interval [r, 



mm 1 ' max 



^pcrihcl 



27r-2 



du 



there is a perihelion shift 



= 27r - Aujik , 



(45) 



where and ei+i are that zeros of P^iu) related to r^ax and r^-iin and the path Ofc surrounds the interval [ei,e,;+i]. 
The actual ujik to be taken also depends on the position of the complex zeros of Pb^u). 



B. Schwarzschild— (anti-)de Sitter in 11 dimensions 



For d ~ 11 and q = we obtain from (1141) the equation of motion 

2 / / 



du \ ' 
dip] 



A 



A 



4 + Au^ - li^ + A(Ai- 1) + — U + — A = 4:Pe{u) 



45 / 45 

and from pH)) the effective potential in the S(a)dS space-time in 11 dimensions 



(46) 



(47) 



In analogy to the case of 9 dimensions we infer from ((^S)) that Pe{u) has two or no real positive roots, thus, leading 
either to escape or to bound terminating orbits. 

For obtaining the analytical solution of (j46p we have to introduce a new variable x by u — ^ + ug , where ug is a zero 

of Peiu). This transforms equation (|46|) containing Pg(u) to a new equation (^x^^ = 4Pi^{x) which is of type (|13ap 

and possesses a solution of the form (|4T|) . The obtained orbits show no qualitative differences to the 9 dimensional 
SdS case, see Figs. [TUl and [TT] for the plots for positive cosmological constant. 



VI. GEODESICS IN REISSNER-NORDSTROM SPACE-TIMES IN 7 DIMENSIONS 

Let us now address the higher dimensional Reissner-Nordstrom (RN) space-time and investigate the effects of the 
charge parameter ij on the timelike geodesies. We restrict ourselves to d = 7 of (fT4| 



d \ 

— I = 4ii (~f]\'^ - Tj^Xu^ + u^ + Xu^-u + A(/i - 1)) = APJu) 
dip J 



(48) 
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(a) pL = 1.008: terminating bound orbit and escape orbit (b) fi = 1.8: terminating bound orbit and escape orbit 





(c) n = 3.13: terminating bound orbit and escape orbit (d) fi = 3.1392073469: terminating bound orbit and escape 

orbit 





(e) /i = 3.13921, terminating 
escape orbit 



(f) /I = 3.5, terminating escape 
orbit 



FIG. 7: Examples of a test particle's motion in 9-dimensional SdS space-time with positive cosmological constant A — 8.7-10 ^ 
and A = 0.2. 



The relation of the zeros of Pe{u) with the parameters A and v are shown in Fig. 12(a) for 77 ~ 0.7. From this plot as 
well as from the effective potential 



1 



1 



1 



(49) 



shown in Fig. |12(b)| it is clear that we have three types of orbits; For one zero (light gray area) we have an escape 
orbit, for two zeros (gray area) we have periodic bound orbit, and for three zeros (dark gray area) we have a periodic 
bound orbit and escape orbit. The relative maximum in Vctf shown in Fig. |12(b)] increases for larger L (smaller A) 
and the relative minimum decreases for smaller q (larger 77). 

In order to obtain explicitly the analytic solution we have to substitute u = — ^ + ue, where ug is the root of 

Peiu), what transforms ((48)) to (^x^^ = 4:Pi^{x). Then the solution of the geodesic equation in Reissner-Nordstrom 
space-time in 7 dimensions is given by the equation 



(50) 
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(a) (A, Ai)-plot, A = -8.7 ■ 10" 




(b) Potential for A = 0.2, 
A = -8.7 ■ 10-5 



FIG. 8: 9-dimensional SadS space-time: (a) (A, -plot (here the dark gray color denotes 3 positive zeros), (b) Effective 
potential. For the corresponding orbits see Fig. [9] 
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(a) fi = 1.008: terminating bound and periodic bound orbits (b) /i = 1.5: terminating bound and periodic bound orbits 



FIG. 9: Orbits of a test particle in 9-dimensional SadS space-time for A = —8.7 • 10 and A — 0.2. 



with the angular coordinate ipz2 from (j33ap and >foc,z2 from p4ap . Xm again is an initial position which can be chosen 
to be a zero of Psix). 

RN space-times have two distinctive features. The first is the existence of two horizons, an inner (Cauchy) and an 
outer horizon which for d = 7 are given by 



1 



^outcr 



-^(2 + 2^13^ 



(51) 



and the second is an antigravitating potential barrier preventing a particle from falling into the singularity which is 
related to lim^^^o ^cff = +oo. For 77=-^= the two horizons coincide and we obtain a degenerate RN space-time. In 



this case the extrema of Vcft in are given by fi = 2 ^Z** and f2 = y A + y + |. At fi the effective potential 

possesses a minimum with Vcsifi) — 0, and at f2 a maximum. That means in particular that a particle with energy 
fi — moves on a circular orbit with radius f = r-mnci — ''outer corresponding to the two coinciding horizons. For 
< /i < 1 a particle necessarily follows a periodic bound orbit. For 1 < fi < T4ff(^2) we have a periodic bound and 
an escape orbit and for fi > Vcff (7^2) only escape orbits. 

For arbitrary rj the extrema of VcS are defined by = — 22:'^A — + Arj'^Xz + 577"* with z = r'^. Using again 
Descartes, there are two or no extrema at positive z. For two extrema there is a range of /i for which there are 
periodic bound as well as escape orbits, and a range of ji with escape orbits only, see the plots in Fig. [131 If there is 
no extremum there are escape orbits only. 

A particle moving along one of the periodic bound orbits shown in Fig. [T3] crosses, say, from the outside, first 
the outer and then the inner horizon, both indicated by black circles. After crossing the inner horizon the particle 
is reflected by the charge induced potential barrier and again crosses the Cauchy horizon in the opposite direction. 
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(a) (A,Ai)-plot, A = 8.7 • 10- 
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M = 4.38 



fj. = 3.0 
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(b) Potential for A = 0.15, 
A = 8.7 ■ 10-5 



orbits see Fig. 1111 
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FIG. 10: 11 dimensional S(a)dS space-time: (a) (X,^) plot (color code as above), (b) Effective potential. For the corresponding 
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(a) /I = 1.0: terminating bound orbit and escape orbit (b) fj, = 4.38: terminating bound orbit and escape orbit 

FIG. 11: Geodesies of a test particle in in 1 1-dimensional S(a)dS space-time with positive cosmological constant A = 8.7- 10"^ 
and A = 0.15. 



thereby entering a new copy of the RN space-time. This can be inferred from the Carter-Penrose diagram of the RN 
space-time shown, e.g., in 41| or [s^. By proceeding further along its f-periodic motion, the particle again and again 
enters new copies of the RN space-time within its analytic continuation. This may be called a many-world periodic 
bound orbit. Since > 0, all bound orbits have to cross both horizons. 

The same happens with the escape orbits shown in Fig. 13(e)|(f) The incoming particle crosses both horizons 
(see the inset in Fig. 13(f) I, is reflected by the potential barrier and, by crossing again the horizons in the transverse 
direction, enters a new copy of the RN space-time. Analogously this may be called a two-world escape orbit. Such a 
scattering or flyby makes the spacecraft disappear into another universe. 

The perihelion shift for periodic bound orbits is 



A. 



pcrihcl 



ei + i 



xdx 



= 27r - 4a;' 



2k 



(52) 



where the zeros and e^+i of P5(x) are again related to the range of motion [r,nin, r-max] and the path ak surrounds 
the interval [ei,ei+i]. Again, the actual 0J2k depends on the position of the other zeros of P5 under consideration. 
Due to the fact that the motion does not take place within one universe it is not clear how one could measure this 
perihelion shift. One possibility might be that an astronaut adds up the angles measured in the various universes 
starting from the position of largest distance with respect to the gravitating body to the next position of largest 
distance. In this case only the astronaut can carry through the measurement and only he or she has access to the 
result. This kind of effect may be called a "many- word perihelion shift". 
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(a) 



(A, M)-plot, »? = 0.7 



(b) Potential for A = 0.35, rj = 0.7 



FIG. 12: 7-dimensional RN space-time: (a) (A,/i)-plot (gray scales as above), (b) Effective potential for r] — 0.7. For the 
corresponding orbits see Fig. 1131 



VII. GEODESICS IN REISSNER-NORDSTROM-(ANTI-)DE SITTER SPACE-TIMES 



As a last example of the motion in higher dimensional space-times, we discuss RN-de Sitter space-times in 4 and 
7 dimensions. This is the most compound case due to a large number of parameters which includes the mass and the 
charge of the gravitating source and the cosmological constant. Owing to its complexity, we start the investigation 
for 4 dimensions and then consider 7 dimensions. 



A. Reissner— Nordstrom— (anti-)de Sitter in 4 dimensions 



Since already the solution of the geodesic equation in 4-dimensional S(a)dS space-times requires the application 
of hyperelliptic integrals [13, |3^ the same is true for the geodesic equation in 4-dimensional RN-(anti-)de Sitter 
space-times 



f-^ + (3X{n - 1) + a) f'' + 3Af^ - 3{f]X + l)f2 + 3r - 3r; j = ^Peir) 



(53) 



The zeros of (1551) are represented in a {X, fi)-plot Fig. 14(a) This is in line with the form of the effective potential 

(54) 



1 ^^2 , V 



shown in Fig. [HT b) for A 
The substitution f = - 
with the general solution given in (E 



0.1, T] ^ 0.1, A = 8.7- IQ-^ 
+ rg, where rg is a root of Pe{f 



reduces ^ to (^xj^^ = which is of type (fT3a|) 

There are periodic bound orbits presented in the first part of Fig 



where particles cross both horizons and then enter another universe 
second part of Fig 



15(a) -(c) 



15(a) -(c) 



Corresponding escape orbits are shown in the 
In the second part of the Fig. 15(a) the influence of the cosmological constant appears 



as a barrier for particles coming from infinity. Further two-world escape orbits for particles coming from infinity and 
crossing the horizons are illustrated in Fig. 15(d)|(e)" for different values of /i. 

A particular property of VcS is shown in Fig. |14(c) for A — 0.3576404 which corresponds to the black area in 



Fig. 14(a; 



The figure reveals three regions where a test particle can move, and in two of them orbits are bound 
are shown in Fig. and exhibit a many-world periodic bound orbit in Fig 



The geodesies for /i = 
ordinary periodic bound orbit in Fig. |16(b)[ and an escape orbit in Fig 



16(c) 



16(a) 



an 



Since in general there are two periodic bound orbits we may obtain two perihelia shifts. Each is given by 



A 



ihcl 



27r- 2 



xdx 



VK{x)/3 



= 2n - iuj' 



2k 



(55) 
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(a) fi = 1.01, r] = 0.2: many— world periodic bound orbit and (b) /i = 1.82, rj = 0.2: many— world periodic bound orbit and 
escape orbit escape orbit 




(c) n = 1.1, r] = 0.7: many-world periodic bound orbit and (d) fj, = 1.6, r] = 0.7: many— world periodic bound orbit and 
escape orbit escape orbit 




(e) ij, = 1.9, rj = 0.2, (f) fj. = 1.9, rj = 0.7: 
two— world escape orbit two— world escape orbit 



FIG. 13: Orbits for a test particle motion in Reissner-Nordstrom space-time in 7 dimensions for A = 0.35. 



where the zeros and ei+i of P5{x) again correspond to the considered range of motion [7')^^;^' ''max] and the path at 
surrounds the interval [ei,ei+i]. In our case j € {1} (only one bound orbit) or j e {1,2} (two bound orbits). While 
in the 4-diinensional RN space-time for /i < 1 the two periods related to two periodic bound orbits are the same, the 
appearance of the cosmological constant removes this degeneracy so that the two orbits possess different periods. 

It is again interesting to investigate naked singularities. In the 4-dimensional RN space-time we have two horizons 
for T] < J (or q < M). For t] — \ both horizons become degenerate and for r] > j a naked singularity appears. 

The horizons for a 4-dimensional RN-de Sitter space-time are given by the real positive zeros of the equation 

^f* - f2 + r - 77 = . (56) 

o 

From Descartes' rule we infer: For negative A there are two or no positive zeros. For positive A there are three or one 
positive zeros, one of these is the cosmological horizon. For both signs of A two zeros may combine to a degenerate 
horizon. 

It is possible to state the algebraic condition on 77 as a function of a given A for which a degenerate horizon occurs. 
In Fig. [17] we show the effective potential ([M)) for two values of the charge parameter, rj = 0.26 and 77 = 0.33; Fig. [18] 
presents the corresponding orbits. The minimum of VoS is negative for small 77, it vanishes for the critical rj, and it lies 
above the f-axis for larger values of ?]. For some larger 77 the minimum disappears. Because of the potential barrier 
a test particle does not fall into the singularity, but in contrast to the many-world periodic bound orbits discussed 
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(a) (A,At)-plot, A = 8.7 • 10-5, ^ = 0.1 
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(b) Potential for A = 0.1 
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(c) Potential for A = 0.3576404 



FIG. 14: 7-dimensional RN-(anti-)de Sitter space-time: (a) (A, /i)-plot (the black region indicates 5 positive roots, dark gray-3 
and light gray-1). |(b)| Effective potential for ri — 0.1 and A — 8.7 ■ 10^^. For the corresponding orbits see Figs. llSl and [161 





(a) /i = 0.95: many-world periodic bound orbit and escape (b) fi = 1.3: many— world periodic bound orbit and escape 

orbit orbit 





(c) /i = 2: many-world periodic bound orbit and escape orbit (d) fi = 2.1: two-world 

escape orbit 



(e) /I = 5.0: two-world 
escape orbit 



FIG. 15: Orbits of a test particle in RN-de Sitter space-time in 4 dimensions for A = 0.1 and rj = 0.1, A = 8.7 • 10 



above the particle will stay in the same universe (compare Fig. 18(a) and e.g. Fig. 16(a) I 



B. Reissner— Nordstrom— de Sitter in 7 dimensions 



For the 7-dimensional RN-de Sitter space-time equation (|14p becomes 



(57) 
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(a) many-world periodic bound (b) periodic bound orbit (c) escape orbit 

orbit 



FIG. 16: Three geodesies in RN-de Sitter space-time in 4 dimensions for the parameters A — 0.3576404, p — 0.888, rj — 0.1, 
and A = 8.7- 10"^ 




FIG. 18: Geodesies in the RN-de Sitter space-time in 4 dimensions with naked singularity for rj = 0.26 and rj = 0.33, 
A = 8.7 • 10"^ 



The number of positive zeros of Pe{u) (or Puir)) is drawn in Fig. 19(a) and the effective potential 

^4 ' 



cff 



1 



1 




(58) 



for A = 0.05, 77 = 0.4 and A = 8.7 -10 in Fig. 19(b) The dark gray area of 3 real positive roots corresponds to 



many-world periodic bound orbits and escape orbits shown in Fig. [201 and the gray area to two-world escape orbits. 
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Introduction of a new variable x such that u ~ + uq, where uq is the root of Pe{u), reduces ([57]) to (^x^^ 
4P5 (x) which is of type (|13ap with the solution given in ([50]) . 




(a) (A,^j)-plot, A = 8.7 ■ 10^5, = 0.4 
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(b) Potential for A = 0.05, 
A = 8.7- 10~5, T] = 0.4 



FIG. 19: 7-dimensional RN-de Sitter space-time: (a) (A, /j) -plot (gray-scales as above), (b) Effective potential for -q = 0.4 
and A = 8.7 • 10~^. For the corresponding orbits see Fig. 1201 
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(a) fi = 1.1: many-world periodic bound orbit and escape (b) fi = 4.0: many— world periodic bound orbit and escape 

orbit orbit 
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(c) fi = 8.4: many-world periodic bound orbit and escape 

orbit 



FIG. 20: Orbits of a test particle in RN-de Sitter space-time in 7 dimensions for A = 0.05 and 77 = 0.4, and A = 8.7 ■ 10 



For a periodic bound orbit in the interval between points and e^+i and the path surrounding this interval the 
perihelion shift Aporihoi is defined as 



A 



perihcl 



27r-2 



xdx 



2tt - 4w2fc 



(59) 
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VIII. LIGHT RAYS 



For light rays e = and the geodesic equation (fT2)) reduces to 



dfV.f 2A / 1 n''-^ , . , 



the structure I ^ ) = Pd-i[ 



dip J \ {d- 2) 

Here the cosmological constant only contributes to the energy parameter, it does not change the mathematical 
structure of the differential equation. The only remaining parameter of relevance is rj. For rj = the equation has 

' — ^- i(r)/f''^'^ and for 77 7^ we have ( — P2{d-2){f) / ^'^'-'^^ ■ For a space-time with odd 
dimension d we get for 

The effective potential for light is 

^cff = f{r)^ = ~ A(d - - 2) ^ Af2 " ^ Af2(rf-2) " 

It is clear that the cosmological constant does not play any role in the motion of light rays. For higher dimensions 
the attractive —l/f'^^^ becomes more pronounced. Also the repulsive part related to the charge becomes stronger for 
higher dimensions and also for larger charges. 

In the case of Schwarzschild and S(a)dS space-times the equations of motion are 



dfV III, 2A 



here are two sign changes in the polynomial on the right-hand-side indicating a terminating bound and an escape 
orbit, or a terminating escape orbit. This can by verified by the form of the effective potential. The equation of 
motion can be analytically solved with elliptic functions for 4, 5, 7 and with hyperelliptic functions for 6, 9 and 11 
dimensions. No analytic solution is known for 8, 10, 12 and all higher dimensions. 
In a RN or RN(a)dS space-time the equation of motion has the form 

f dfV 1 ( ( , 2A " 



- 



+ f2(d-2) _ ~2(rf-3) ^ ~d-3 _ (g4) 



\dip J f2(d-i) n {d-i){d 

If the first term is positive then Descartes' rule states that there are three or one real positive zeros of the polynomial 
implying the existence of many-world periodic bound and escape orbits or two-world escape orbits. That means light 
can disappear into another universe or appear from another universe. For a large negative A escape orbits no longer 
exist. This can also be seen from the corresponding effective potential. Using (|6ip the equation of motion can be 
solved for 5 and 7 dimensions, in 5 dimensions by elliptic functions, in 7 dimensions with hyperelliptic functions as 
shown above. 

IX. RADIAL MOTION 

For radial motion, i.e., for L = 0, the geodesic equation ^ for point particles reduces to 

drV _ 1 2Af2 

Ts) - ^ " ^ + ^ + (d - l)(d - 2) " ■ (65) 

This equation can be integrated along the lines presented above for non-radial motion. The effective potential for 
this radial motion is 

1 2Af2 Tyrf-s 

^eff - 1 - ^ - [d-l){d-2) ^ ^^^^ 
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Of interest are the equilibrium positions given by the vanishing of the first derivative of the effective potential 

= (d _ 3^-^ _ __i^__,~2(.-2) _ 2(rf _ 3)^.-3 (67) 

If 7] = then there is a solution only for A > which, however, corresponds to a maximum and, thus, gives no stable 
position. This equation has stable solutions only ii rj ^ 0, independent of the dimension. For simplicity, we now 
choose A = 0. Then (|67l) has the solution fo — with the effective potential Vcs_{fo) = 1 ^ \'n^~'''- Therefore, in 

any dimension for -q > a charged solution allows particles to stay at rest at a stable position at a finite value of 
the radial coordinate. 



X. CONCLUSION AND OUTLOOK 



A given gravitational field can only be analyzed and interpreted through the exploration of the geodesies of particles 
and light rays. In this paper we have discussed the motion of test particles and light rays in higher dimensional 
space-times of spherically symmetric gravitational sources endowed with mass and electric charge and a cosmological 
constant. We discussed the general structure of the resulting orbits. One result of this analysis is that only in 4- 
dimensional Schwarzschild and Schwarzschild~de Sitter space-times it is possible to have stable periodic orbits which 
do not cross horizons. In Reissner-Nordstrom space-times of any dimension we have stable orbits which, however, 
periodically disappear into other universes. This raises questions about the operational realization and interpretation 
of notions like the perihelion shift and scattering angles. Furthermore, due to the structure of the angular momentum 
barrier and the cosmological force only in 4 and 5 dimension we have a rich variety of orbits. For 6 and higher 
dimensions we always encounter the same two types of orbits. 

We analytically integrated the equations of motion for test particles and light rays in Schwarzschild and 
Schwarzschild-de Sitter space-times in 9 and 11 dimensions, in Reissner-Nordstrom space-times in 7 and Reissner- 
Nordstrom-de Sitter space-times in 4 and 7 dimensions by applying a method based on a solution of Jacobi inversion 
problem restricted to the set of zeros of a theta function and described in 34, 35]. The explicit integration was 
possible because in these cases the underlying polynomial appearing in the equation of motion is at most of order 6. 

The next step would be to find the solutions of geodesic equations with an underlying polynomial of 7th and higher 
order (cases indicated with "-" in Table [l|. In these cases we have to enlarge the number of variables to be three or 
more. It is not clear currently how to constrain the Abel mapping between 3 or higher dimensional spaces in order to 
reduce the number of variables. 

Meanwhile, it is of special importance to find analytical solutions of the equations of motion. Because of the, in 
principle, arbitrary high accuracy of analytical solutions one can test with these solutions numerical codes, e.g., for 
the dynamics of binary systems in EMRIs (Extreme-Mass-Ratio Inspirals). Another application is that the offered 
solution can serve as starting point for an advanced approximation scheme in more complicated situations appearing 
in stellar, planetary, comet, asteroid or satellite dynamics. 
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